Cavity state preparation using adiabatic transfer 
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We show how to prepare a variety of cavity field states for múltiple cavities. The state preparation 
technique used is related to the method of stimulated adiabatic Raman passage or STIRAP. The 
cavity modes are coupled by àtoms, making it possible to transfer an arbitrary cavity field state from 
one cavity to another, and also to prepare non-trivial cavity field states. In particular, we show how 
to prepare entangled states of two or more cavities, such as an EPR state and a W state, as well as 
various entangled superpositions of coherent states in different cavities, including Schrodinger cat 
states. The theoretical considerations are supported by numerical simulations. 
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I. INTRODUCTION 

A recent paper pj presented an efficient method to 
adiabatically transfer field states between two different 
cavities. The scheme is closely related to stimulated Ra- 
man adiabatic passage, or shortly STIRAP STI- 
RAP was first used to coherently control dynamical pro- 
cesses in àtoms and molècules. Two external làser pulses 
drive population between an initial and a final state in 
an atom or molecule, through an intermediate level. One 
pulse couples the initial state to the intermediate state 
and the other pulse couples the intermediate and final 
state. The pulses are applied in a counterintuitive way, 
in the sense that the pulse that couples the final and in- 
termediate states is turned on first. The pulses do have 
to overlap though, and in order for the process to work 
successfully it has to be adiabatic, as the name suggests. 
Population will then follow the instantaneous eigenstates 
adiabatically. One of the eigenstates is of particular in- 
terest, namely the dark state. This state has eigenvalue 
zero, and the intermediate state is never populated dur- 
ing the evolution. 

In the method suggested in 0, a two- level atom in- 
teracts with two cavities. In this scheme, the couplings 
between the atom and the two cavities correspond to the 
the two làser pulses in traditional STIRAP. As the atom 
traverses the cavities it will see the varying shape of the 
mode it interacts with, and consequently, the coupling 
becomes time-dependent. By letting the cavities partly 
overlap spatially, it is possible to realize a situation very 
similar to STIRAP. In fact, if the state, adiabatically 
transmitted between the cavities, is a one photon state 
|1), the corresponding Hamiltonian (in the dipole and 
rotating wave approximations) looks exactly the same as 
the Standard STIRAP one. The ingenious feature of the 
method is that it works for any field state, not just the 
one photon state. The Hilbert space will, of course, in- 
crease when larger photon number states are involved, 
and therefore the adiabaticity constraints become more 
stringent There is still a dark state with zero pop- 
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ulation in the upper atòmic level, even for general field 
states. 

Other schemes, where the atom experiences a varying 
mode shape as it traverses the cavity, have also been sug- 
jested_for_ adiabatic state preparation of the field modes 
:ver, these 
example, ii 

type atom is used, in [íj Q, H, a strong external classical 
làser field is utilized and in Q only one cavity and one 
two-level atom is considcred. 

In this paper we will extend the model in p| to more 
complex systems involving more than just one two-level 
atom and two cavities. As we have mentioned, in the 
one photon case the model in is analogous with the 
traditional STIRAP. Likewise, the extensions made in 
this paper are related to similar generalizations of the 
traditional STIRAP, if we consider the one photon case. 
General situations for multi-level STIRAP has been an- 
alyzed in several papers; just to mention a few, see 
ES IH El El El ■ By including more àtoms and cav- 
ities, we will show that various interesting field states 
can be prepared. Duc to the fact that the dimension 
of the accessible Hilbert space easily blows up when the 
photon number is increased in these extended models, 
we will choose the transferred field state to contain just 
one photon in our numerical simulations. However, in 
the adiabatic limit, the system is solvable also for higher 
photon numbers. Using more photons only means that 
the adiabaticity constraints are stricter, as mentioned 
above. As compared with the method in El> we will 
note that also these more complicated systems have an 
adiabatic dark state, which will be used for the evolution. 
It will be shown that it is possible to entangle spatially 
separated cavities, and prepare, for example, EPR or 
W field states, but also more complex entangled states. 
By making atòmic measurements, it is feasible to create 
Schrodinger cat states. The setups given in this paper 
are only a couple of examples, and others are of course 
possible; we just illustrate the bàsic idea. We consider 
preparation of the various field states, but the methods 
could equally well be applied for creating different atòmic 
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states if desired. 

The outline of the paper is as follows: In section [Q] 
we review the bàsic idea and properties of the method 
presented in We introduce the adiabatic eigenstates 
and explain the dynamics behind the transfer of arbitrary 
field states between two cavities. In section lïïïl we con- 
sider two different setups, which we call the "H" configu- 
ration, consisting of three cavities and the "star" config- 
uration, which could contain any number M of cavities. 
In the H configuration we show how a state is transfered 
between two spatially separated cavities by virtual pass 
through a third cavity and it is also explained how EPR 
states could be prepared. The other model, the star con- 
figuration, could also be used for achieving EPR states 
as well as W states and generalizations of these states. In 
section lïvl we make use of a third atòmic level and pro- 
jective atòmic measurements for preparing various types 
of Schròdinger cat states. Finally we conclude with a 
summary and discussion in section Ivl 



II. ADIABATIC TRANSFER BETWEEN 
CAVITY MODES 

We will first briefly review how to adiabatically trans- 
fer a quantum state from one cavity mode to another, fol- 
lowing |íj. We consider a situation where there are two 
cavity modes interacting with a single two-level atom. 
The Hamiltonian for this system is a generalisation of 
the widely used Jaynes-Cummings model 15], 

H = —Lo{a z + 1) + í\à{ài + f^aíjé^ 

+ Oiói + 92à 2 )a+ + Í9ià\ +g 2 àl)a~. (1) 

Here à\ and à\ are the boson creation operators for cav- 
ity modes 1 and 2, respectively, <7 Zl a + and er~ are the 
Pauli z and the raising and lowering operators for the 
atom, and gi(t) and g 2 (í) describe the time-dependent 
coupling between the light and the two-level atom. The 
basis states for the system are of the form 



\ni,n 2 , 



\ni)\n 2 )\s) 



(2) 



where n\ and n 2 refer to the number of excitations in 
mode 1 and 2, and s = ± refers to the state of the 
two-level atom, with a z \s) — s\s). In the following 
we will assume that the cavity modes are degenerate, 
íii = Cl 2 = íl, so that perfect transfer of excitations be- 
tween the modes is possible. If we start with a single 
excitation in mode 1 and the atom in its ground state, 
then the accessible Hilbert space is spanned by the three 
states 



|1,0,-),|0,0,+),|0,1,-). 
The Hamiltonian commutes with the operator 



N 



1 



(cr z + 1) + à\ài + à\à 2 , 



(3) 



(4) 



so that we can work in an interaction picture, with the 
Hamiltonian 



H' 



H-CIN 
A{a z + 1) 



(5) 



[(gi(í)ai + g 2 {t)à 2 )<7 + + h.c. 



where A = (oj — f2)/2. The atom does not need to be on 
resonance with the cavity modes, i.e. A can be nonzero. 

As in the case of adiabatic transfer between atòmic 
states |2J, , there is an eigenstate of this Hamiltonian 
with eigenvalue zero, given by 



K12 [ff 2 (í)|l,0, 



5i(í)|0,l, 



(6) 



where the normalisation constant is given by 
g\(t) + g 2 (t). Consider the case when 



lira 



9i(t) 
92 (t) 

hm ^ 



(7) 



If the couplings g\ (t) and g 2 (t) change slowly enough, the 
system will start in the state |1, 0, — ), and end up in the 
state |0, 1, — ), following the adiabatic eigenstate given in 
equation ijjjj. This method is called stimulated Raman 
adiabatic passage or STIRAP 0, Q . The exact shapes of 
the pulses g\ (t) and g 2 (t) do not matter, as long as they 
vary slowly enough and conditions @ hold. The pulse 
sequence is counterintuitive in the sense that the two ini- 
tially empty lcvels are coupled first, and only then is the 
initially populated level coupled to the "middle" level. 
The two pulses gi(t) and g 2 (í) must, however, overlap. 

By choosing lim t _ 00 g 2 (t)/gi(t) = 1 instead of 0, we 
can also adiabatically reach the state 



-L(|i,o,-)-|o,i,-)), 



(8) 



or, by choosing another suitable ratio between gi(t —* 00) 
and g 2 (t — ■+ 00), we can reach any superposition of 
|1, 0, — ) and |0, L — ). This process is referred to as frac- 
tional STIRAP Q. 



A. Transfer of an arbitrary cavity field state 

Also more than one field excitation can be transferred 
between the cavity modes pj. For example, a Fock state 
|n) in mode 1 can be transferred to mode 2. We can write 
the adiabatic state © as 

|* od )=it| ,0,-), (9) 

where the boson operator A> is defined as 

i 1 =K 12 (g 2 à\ -gia\). (10) 

The Hamiltonian, on the other hand, can be written as 

H' = A{a z + l)+K^ 1 (Ba+ + B^a-), (11) 
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where the boson operator B< is given by 

Èï =K 12 ( gi à\+g 2 àt). (12) 
We find that [È, À^] = 0, so that the state 



i^) = W 7i(^ t mo.-) 



(13) 



is an adiabatic state, since H'\^" d ) = 0. Choosing the 
couplings so that conditions J7J| hold, we immediately 
find that the state |n,0, — ) adiabatically changes into 
|0,n,-). 

More generally, we can consider the adiabatic state 

/(At)|0,0,-)=C n ^-^|0,0,-). (14) 

If the couplings again satisfy conditions Q, and if we 
choose the pulses so that gi/g 2 < 0, then the state 
f(à\)\0, 0, — ) will adiabatically change into f(àl)\0, 0, — ). 
For example, a coherent state \a) can be transferred from 
cavity mode 1 to cavity mode 2 by choosing 

\*ad) = cxp (-^) exp(ait)| 0j 0, -). (15) 



III. ADIABATIC TRANSFER WITH MÚLTIPLE 
CAVITIES 

A. Three cavities and two àtoms in an "H" 
configuration 




cavity 1 



FIG. 1: A possible setup of the three cavities (1, 2 and 3) 
and the two atòmic (o and 6) trajectòries for the "H configu- 
ration" . 

We will now move on to consider cavity state transfer 
in a situation where we have three cavities and two àtoms. 
Suppose cavities 1, 2 and 3 are placed so that cavity 1 is 
overlapping with both cavities 2 and 3. Atom a is placed 
in the crossing between cavities 1 and 2, and atom b in 
the crossing between cavities 2 and 3, as shown in figure 
^ The Hamiltonian for this system is given by 



H 



\bJ b {<J bz 



+ il 2 à [ 2 à 2 + ^ 3 4 ó 3 + [(giaài + g 2a à 2 )<r+ (16) 
+ (duài + gzbà?,)^ + h.c] , 



where cr a (f))zi a t(b) an< ^ °a(6) re f er to atom a(b), and à\ 
and ài are the creation and annihilation operators for 
cavity i. We have denoted the coupling strengths between 
cavity i and atom a as gi a , and correspondingly for atom 
b. The number of excitations in the systems is conserved, 
and we find that the Hamiltonian commutes with the 
operator 



N=\{a az 



1/ 

2 (CTq 



D 



à\à 2 



(17) 



In the following we will assume that Í2i = Cl 2 = ÍI3 = Cl. 
Otherwise perfect transfer of cavity field states would not 
be possible, since energy is conserved. In the interaction 
picture, we form the Hamiltonian 



H = H-ÜN = A a (cr az + 1) + A b (a bz + 1) 



[(giaài 



92aà 2 )a+ 



(guài +93bà 3 )a^ 



h.c. 



(18) 

where A Q = (ui a — fl a )/2, and similarly for b. We now 
write the basis states as \rii, n 2 , 713, ± a ,±f,), where the 
three first entries refer to the number of photons in 
cavities 1, 2 and 3, and the two last entries to the states 
of the àtoms. The subspaces with exactly one excita- 
tion in the system is spanned by the five basis states 
|0, 1,0,-,-), |0, 0, 0, +, -), |1, 0, 0, -, -), |0, 0, 0, -, +) 
and 1 0, 0, 1 , — , — ) . Using this ordering of the basis 
states, the Hamiltonian in matrix form for this subspace 
becomes 



H 



( g 2a \ 

g* 2a A a g la 1 

g* la g lb 

g* b A b g 3b 

V g* b / 



(19) 



This Hamiltonian has an adiabatic eigenstate with eigen- 
value zero. Making the Ansatz (C 2 , 0, Ci, 0, C 3 ) T for 
this state, the condition on the coefircients Ci becomes 
9 2a ^2 + giaCi = g{ b C\ +g3 b C b = 0, so that the adiabatic 
eigenstate is 



K{gi a g3b, 0, -g 2a g 3 b, 0, g^g^) 1 



(20) 



where K is a normalisation constant. We see that there 
should be a possibility of transferring the state of cavity 2 
directly to cavity 3 with very little population in cavity 1. 
For a thorough exposition of adiabatic transfer between 
atòmic levels with múltiple intermediate states, see |l0j . 
The theory can be directly applied to cavity state transfer 
as well. To achieve transfer from cavity 2 to cavity 3, we 
should start with 



and finish with 



|5laff36| > \gib92a\, 



lsi&52a| > \gia9ab\, 



keeping 



\giag3b\ 2 + \gibg2a\ 2 > |52a33b| S 



(21) 
(22) 
(23) 
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all the time. There are many possible pulse sequences 
satisfying these conditions. A few possible coupling se- 
quences will be discussed in the next subsection. In all 
cases we start with one field excitation in cavity 2. 

As for the case where two cavity modes are coupled by 
one atom the transfer of arbitrary cavity states from 
mode 2 to mode 3 will also be possible. If we form the 
"adiabatic operator" 

À\t) =K{t)[g ía {t)g 3h {t)àl-gÍ a {t)g 3h {t)à\ + gl b {t)g* 2a {t)àt\ , 

(24) 

where K(t) is a normalisation constant, then, in the 
adiabatic limit, if we start in the state /[A^(0)]|0), we 
will also stay in the state /L4^(í)]|0) as the couplings 
are changed. For example, starting in /(à|)|0, 0, — ), we 
can adiabatically transfer the cavity state to mode 3, 
/(0.3) |0, 0, — ) . As before, this means that we can trans- 
fer not only one field excitation, but also, for example, 
number states, where f{A^) = A* n , and coherent states, 
where /(At) = exp (|a| 2 /2) exp(ait). 

B. Numerical simulations of the "H" configuration 

For all the numerical simulations in the paper we use 
Gaussian pulses for the couplings, of the form 

9iAt) = Giue X p(- ^ t ~l iv)2 ). (25) 

The index i stands for the í'th cavity and v for atom 
V) cavities will be labeled with numbers and àtoms with 
letters. If there is only one atom present the atòmic index 
will be omitted. G is the coupling amplitude, and it 
will be chosen the same for all pulses in the different 
examples, except for a couple of examples in the next 
section. The indices will be omitted when the G:s are all 
the same. The parameter íj„ gives the pulse center and 
the width is given by <7j„ . We are using scaled parameters 
with h = 1. Time t and the pulse widths a are given in 
units of a suitable characteristic time T, and G and A in 
units of KT~ X . 

We will consider two possible pulse sequences for adia- 
batic transfer in the "H" configuration. The first pulse se- 
quence, which is shown in figurc[21 is completely counter- 
intuitive, in the sense that we start by coupling cavity 3 
and atom b, then cavity 1 and atom b, followed by cav- 
ity 1 and atom et, and finally cavity 2 and atom a. This 
could for example be achieved if the cavities are crossing 
each other horizontally, partly overlapping, and we let 
the atom b traverse first cavity 3 and then cavity 1, and 
similar ly for atom a and cavities 1 and 2. The parameters 
in the figure are t^b — —5.22, tu, = —1.72, t\ a — 1.78 and 
t 2a = 5.28, cr = 3, A = and G = 100. The dynamics is, 
for A = 0, determined by the dimensionless adiabaticity 
parameter Ga 0. 

The pulses are seen in the left plot and the popula- 
tions in the right one. As shown in figure [21 numerical 




-20 -10 10 20 -20 -10 10 20 



t t 

FIG. 2: The figure to the left shows our first example of 
a pulse sequence for realizing complete population transfer 
from cavity 2 to cavity 3 with minimal population in the in- 
termediate cavity 1 for the "H" configuration. The pulses 
are ordered in a completely counterintuitive way, from left to 
right 336, (/it, gia and gi a - Time is given in units of a suitable 
characteristic time T . The widths of the pulses are all a = 3, 
also in units of T, and the maximum amplitudes are G = 100 
in units of T -1 . The other plot shows the populations \Ci\ 2 
(i = 1, 2, 3, 4 or 5) as a function of the scaled interaction time 
t. It is clear that population is transfered adiabatically from 
the second cavity (sòlid line marked |Ci | 2 ) to the third cavity 
(dotted line marked |Cs| 2 ), without remarkable population in 
cavity 1. The final population in the third cavity is 99.8 %, 
and maximum population of cavity 1 during the process is 0.2 
%. 



simulations confirm that an excitation in cavity 2 can 
be transferred adiabatically to cavity 3, while the pop- 
ulation in cavity 1 remains small in between. The final 
population in state |0, 0, 1, — , — ) is 99.8 % and maximum 
population in cavity 1 is 0.2 % and is located around 
t = 0. The coupling amplitudes are rather large in this 
example in order to have an adiabatic process and corre- 
spondingly a successful transfer. This is due to the fact 
that the population virtually passes through three lev- 
els, |1, 0, 0, — , — >, |0,0,0, +,-) and |0,0, 0, -, +), instead 
of just one in the Standard STIRAP. However, it is still 
clear that if the procedure is slow enough it is possible to 
transfer the population adiabatically. It is also possible 
to switch the order of the two middle pulses ^3 • 

In this example, the population transfer takes place 
mainly when all four pulses differ from zero, when the 
product g prod = gi a g2agib9sb ^ 0. Letting g prod increase 
by making the pulses overlap more in time, it is possi- 
ble to have efficient population transfer from state one 
to state five with a smaller adiabaticity parameter Ga. 
However, the price one has to pay is that in this case, 
the intermediate states become more populated during 
the evolution, since condition i|23|) is not as well satis- 
fied. Thus, there is a tradeoff between strict adiabaticity 
constraints (large Ga) and small population of interme- 
diate states, or weaker constraints but population of the 
intermediate states during the transfer. 

Another possible coupling sequence is shown in figure 
Here the coupling time between atom a and cavity 1 
is longer than the coupling time between and atom b and 
cavity 3, and these two couplings are centered around 
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FIG. 3: The same model as in figure |3] but with the second 
choice of pulse sequence, where the pulses are allowed to have 
different widths. The pulses, shown to the left, come in the 
following order: first g\ a (sòlid) and g^b (dotted) at ti a = 
Í3(, = —3 and then g lb (sòlid) and g^ a (dotted) at tu = Í2a = 
3. The widths for gi a and gu (sòlid curves) are a = 6 and 
for the other two pulses (dotted) a = 2, and the maximum 
amplitudes are as in the previous example G = 100. Time 
is given in units of T and the pulse height in units of T~ . 
The population transfer, shown in the right plot, is similar 
to the previous example, with a final population in cavity 



3 ICsl = 99.; 



and a maximum population in the middle 



cavity equal to 0.8 %. 



the same time. Similarly, the coupling time between 
atom ò and cavity 1 is longer than the coupling time 
between atom a and cavity 2, and these couplings are 
also centered around the same time. This coupling se- 
quence also satisfies the conditions (|21I - I23[1 . It could be 
achieved by making the diameters of làser beams 2 and 
3 smaller than the diameter of làser beam 1 . Numerical 
simulations confirm that this coupling sequence works, 
and the population is transferred from cavity 2 to cavity 
3 with very little population of cavity 1 during the trans- 
fer. The parameters in this second choice for the pulses 
are, t\ a = Í3& = —3 and tu, = t% a = 3, <J\ a = <J\b = 6 and 
(T2a = o"3h = 3 and again G = 100. The plot to the right, 
for the population transfer, looks similar to populations 
in figure [21 and here we have final transfer in cavity 3 
| C5 1 2 = 99.8 % and maximum population in the middle 
cavity 1 | C3 1 2 = 0.8 %, thus a small fraction more than 
in the previous example. 

Since cavity 1 remains almost unpopulated for the cou- 
pling sequences we have discussed, relatively large losses 
in cavity 1 should not affect the emeieney of the state 
transfer. This is also confirmed by numerical simula- 
tions. In order to investigate the effect of losses in the 
intermediate cavity we add a loss term 

6 = e'^ (26) 

to the derivative of the amplitude of the state 
|1, 0, 0, — , — ). To check the advantage of our model, with- 
out population in cavity 1 , compared to a situation with 
population in cavity 1, we simulate a situation were atom 
a transfers the photon first to cavity 1 from cavity 2 
and then atom b takes it to cavity 3. This amounts to 
two consecutive ordinary STIRAPs with population in 
the middle cavity. First we show the population trans- 
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FIG. 4: This figure shows the effect of losses in the interme- 
diate cavity 1. The setup is as a double STIRAP, the first 
atom transfers the photon from cavity 2 to cavity 1 and íi- 
nally the second atom brings it into cavity 3. Note that the 
pulses of the two STIRAP overlap, thus the middle cavity is 
never fully populated. The pulse parameters are G — 100 



(in units of T 



tia 



"3, Í2a 



-1, í 3i) = 1 and tib 



and cria = (72a = <J lb = a 3b = 2 (in units of T). The left 
plot shows the populations without losses in cavity 1, while 
in the right figure, cavity 1 has a decay rate 7 = 0.1. The 
final population transfer from cavity 2 to 3 is reduced from 
100 % to 20 %. This should be compared to, for example, 
using the pulse sequence of figure [5] with losses. If we add 
the same decay rate 7 = 0.1 for cavity 1 in that process, the 
population transfer goes down from 99.8 to 99.7 %. 



fer without losses in cavity 1 in the left plot of figure 01 
and then we add the loss term l|26l) to the Hamiltonian 
with a decay rate 7 = 0.1 and we see the result in the 
plot to the right, the transfer efficiency goes down from 
100 % to 20 %! When adding the same loss term to the 
example in figure [21 the decrease in population transfer 
is only 0.1 percentage units. The parameters for figureQ] 
are t\ a = —3, Í2a = —1, Í36 = P an d t%b = 3, a = 2 and 
G = 100. If we increase the decay rate to 7 = 1, keep- 
ing all other parameters the same, the population goes 
down to 99.0 % in our first method, while in the second 
model, when cavity 1 is populated, no population ends 
up in cavity 3. 

Losses will, however, broaden the lineshape of the cav- 
ity. If the cavity is too long, factors exp(ifcr) coming from 
the propagation in the cavity will most probably disturb 
the adiabatic transfer process, since the line is broadened 
and therefore not only one value, but vàlues of k in an 
interval are involved. For a long lossy cavity 1, the ef- 
nciency of adiabatic transfer from cavity 2 to cavity 3, 
trying to avoid the lossy cavity 1, will be lowered. 



C. Preparation of an EPR state in the "H" 
configuration 

So far we have only been discussing transfer of a field 
state between two cavities, separated in space, but the 
model could also be used for creating entanglement be- 
tween the cavities. Here we give an example of that, and 
the following two sections will consider entanglement in 
more detail. We introduced the adiabatic eigenstate l|20|) 
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with eigenvalue zero, and by choosing the pulses gi v care- 
fully we could transfer population, but, of course, there 
are numerous other interesting pulse sequences. Assume 
that gi a and g^b are turned on simultaneously and then 
gi a and gu are turned on simultaneously. The adiabatic 
state then begins, at t = — oo, as (0, 0, 1, 0, 0) T and ends 
as (— 1, 0, 0, 0, — l) T /v / 2· Thus, by lctting atom a and b 
interact simultaneously with cavity 2 and 3 respectively, 
and then simultaneously interact with cavity 1, the ini- 
tial photon in cavity one will be transfered into an EPR 
statc, 

\EPR) ± = -L(\0,1)±\1,0)), (27) 

of cavities 2 and 3. This procedure is shown in figure 
The pulses are given in the left plot and the populations 
in the right plot. The parameters are ti a — t^b = — 2 
and t ía — t\b = 2, u = 3 and G = 5. Note that here 
the coupling amplitudes (and correspondingly the degree 
of adiabaticity) does not need to be as large as in the 
examples of adiabatic transfer. The photon clearly ends 
up in cavity 2 and 3. That the state is really the pure 
state H27f) . and n °t a mixture, is checked by calculating 
the fidelity between the final state from the numerical 
simulation and the EPR state 

F = \ + (EPR\tp(t = +oo))|. (28) 

With the state obtained numerically with the parameters 
in figure[5l the fidelity becomes F = 0.9999. By control- 
ling the phases of the couplings it would be possible to 
obtain different EPR states. Starting with a general field 
state in cavity 1, the final state would be a more com- 
plicated entangled state of cavity 2 and 3, obtained with 
the method explained in the previous section, by acting 
with the adiabatic operator f{A') on the vacuum. The 
situation is analogous to when a coherent state is split 
by a 50/50 beam splitter. 




t (a.u.) t (a.u.) 



FIG. 5: In this figure it is shown how well the method works 
for preparation of EPR states between cavity 2 and 3. To the 
left we show the pulses, with the parameters G — 5 (in units 
of T" 1 ), tia = 2, Í2a = —2, t 3b = —2 and t lb = 2and a\ a = 
(72a = 016 = °36 = 3 (in units of T). The right plot gives the 
populations, and it is clear that population initially in cavity 
1 (sòlid line) is transfered equally to cavity 2 and 3 (dotted 
and dashed line). Note that in this situation the amplitude 
G is much smaller than in figures [2] and El For the fidelity in 
this example we have F = \(EPR\tp{t = +oo)}| = 0.9999. 




FIG. 6: This figure shows a possible setup for the 'star' con- 
figuration with three cavities. Note that two of the cavities 
should be in the same plane, while one (the initially populated 
cavity) is slightly off the plane. The atom passes through the 
cavities in the middle point of the 'star'. 



D. "Star" configuration 



We can easily extend the situation to more than three 
cavities, or to other setups, such as a ring configuration, 
where the three cavities form a triangle, overlapping each 
others at the còrners of the triangle. In this section we 
investigate a situation with M cavities and one single 
atom coupled to all of the cavities, as shown in figure [ü] 
We will also discuss the effect of adding further àtoms 
coupled to some, but not all, of the cavities. If the atom 
travels along, say, the z-axis, the cavities form a "star" 
in the xy-pl&ne. We assume that M — 1 of them are in 
the same plane, centered around z — 0, and cavity M 
is slightly shifted from z = 0. Initially only cavity M is 
populated and again we take all ÍVs to be identical. 

The effective Hamiltonian for the system is, in the ro- 



tating wave and dipole approximation, given by 



H = A(a a + 1) 



M-l 

gMaàl·l^a + 9a ^ à iCr+ + h.C 
i=l 



(29) 

Note that we have assumed that the couplings are iden- 
tical for the first M — 1 cavities, gi a — g a for i = 
1,2,. ..,M — 1. For simplicity, we again considcr only 
the case with one excitation, N = 1. By labeling the 
states as |1,0, ..,0, — ), |0, 1, 0, — ),..., |0, 0, 1, — ) and 
|0, 0, 0, +), we find the adiabatic eigenstate 

|*)ad = K(-g Ma , -gua, ■ ■ ■ ,-gMa,g a ,0) T (30) 
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t (a.u.) G ( a - u ) A (a.u.) 



FIG. 7: This shows the numerical simulation of the 'star' 
configuration and the preparation of a W state. The left 
plot gives the pulses in time. The parameters are G = 5, 
ai = <T2 = <T3 = <T4 = 2 and ti — Í2 = Í3 = — 1 (dotted 
lines), Í4 = 1 (sòlid line). Time and the pulse widths <r are 
given in units of T and the pulse heights are given in units of 
T _1 . To the right we see the population, and it is easily seen 
that the initial population in cavity 4 (sòlid line) is equally 
transfered to cavities 1, 2 and 3 (dotted lines). The fidelity is 
F = \{W\i/>(t = +oo))| = 99.8%. 

with eigenvalue zero. Thus, if we have 

limt^-oo (s^a) = 0, and lim t ^ +00 (^) = 0, 

(31) 

the photon will be adiabatically transfered from cavity M 
into all other cavities with equal probability and phase. 
With M = 3, the final state in the first two cavities will 
be an EPR state, and with M — 4, we get a so called W 
state, 

\W) =^(|1, 0,0) + |0,1,0) + |0,0,1)). (32) 

For M > 4, it is possible to prepare the natural gener- 
alization of the W state to higher dimensions. A similar 
setup and the generation of W states were discussed in 

m 

In figure 0we show the pulses and populations during 
the passage of the atom, with four cavities (M = 4). 
The parameters are A = 0, G — 5, íi,2,3 = 1, Í4 = 
— 1 and 0+2,3,4 = 2. The dotted lines shows the pulses 
g a (t) and the sòlid line the pulse g4 a {t)- The process is 
counterintuitive like the original STIRAP. In fact, this is 
an "ordinary" STIRAP, but with N— 1 final states, rather 
than just a single one. We clearly see that the population 
is equally split between the first three cavities, and with 
these parameters the fidelity is F = \(W\tp(t — +oo))| = 
99.8. Note that, as for the generation of the EPR state in 
figure[3 the amplitude G is rather small in this example, 
compared to the case of population transfer between the 
cavities in the " H" configuration which is shown in figures 
Hand[2 

Next we show how well the process works for differ- 
ent parameters, changing the coupling amplitude G and 
the detuning A between the atòmic transition frequency 
w a and the common field frequency íí. In figure |Hl the 
parameter dependence of the fidelity F = \(W\ip(t = 



FIG. 8: This figure shows the fidelity F = \(W\ip(t = +oo))| 
as a function of the coupling-amplitude G (left plot) and as a 
function of the detuning A (right plot), for the example given 
in figure|7| In the first plot A = and in the second G — 5 (in 
units of T -1 ), otherwise the parameters are as in the previous 
figure [7| 

+oo))|, in the previous example, is shown; first as func- 
tion of the amplitude G, with A = 0, and then as function 
of A, with G = 5. The other parameters are as in figure 
[7| The fidelity, as expected, increases with the coupling 
and decreases with the detuning. Similar plots could be 
made for the other examples, and the information ob- 
tained would be similar. 

In section [H] we explained how a general Fock state 
\n) is adiabatically transfered between two cavities. The 
same procedure can, of course, be used also in this con- 
figuration. In a similar fashion as in equation (jüj), we 
introduce an "adiabatic operator". Using the pulse se- 
quence above, the adiabatic state l|3U|l will then evolve 
according to 



(33) 

where |0, — ) on the right hand side means vacuum plus 
ground state atom and Ï/N is a normalization constant. 
Here we have also used the multinomial theorem. Know- 
ing how a Fock state transforms, it is easy to calculate 
how a general state in cavity M evolves. States of sim- 
ilar forms as the one above, but for two modes, have 
been discussed for example in 0, . By selecting the 
cocfficicnts in equation (|29() to differ between the indi- 
vidual modes, more general final states can be prepared 
adiabatically. 

In the adiabatic limit the system evolves according to 
the adiabatic states, and the process is robust against 
small changes in the parameters pfij , which is a great ad- 
vantage for example in quantum computing 21, 22J. The 
adiabatic states are, however, sensitive to small changes 
in the Hamiltonian, which will be shown next. If a second 
atom b, also in its ground state, is coupled to only cavity 
j in the " star" configuration during the whole passage of 
the first atom through the M cavities, we have to add an 
interaction term to the Hamiltonian (|29|l of the form 

V = gjbàj&b + h.c. (34) 

We assume that the detuning between the j'th cavity 



8 



and the second atom b is zero, so in the interaction pic- 
ture the atòmic energy vanishes. The shape of g^ is not 
so important as long as it is non-zero during the pro- 
cess. We take it to be constant, but it could also be a 
very broad Gaussian, so that it extends outside the other 
Gaussian pulses g a and gMa, which could be the situa- 
tion if the second atom moves much slower than the first 
atom and only through the j'th cavity. By adding the 
tcrm <|34[) to the original Hamiltonian, the Hilbert space 
dimension obviously increases by one unit, due to the 
state |0, 0, .., 0, — , +), and the corresponding adiabatic 
state i ffity becomes 

\^)ad = K(-g Ma , -gMa, • • ■ ,0, . . . ,-gMa, 9a, 0,0) T , 

where the new is on the j'th position. The added atom 
thus takes away the population in the j'th cavity. In the 
adiabatic limit, the magnitude of gjb is not important, 
just that it is non-zero. In other words, coupling one of 
the 'bare' states in the Hamiltonian weakly to a 'new' 
state drastically affects the adiabatic evolution. If a new 
atom c or atom b is coupled to yet another cavity l during 
the whole interaction, the population of that state would 
become zero. 

The modification in the evolution is shown in figure 
|ÜJ We use exactly the same example and parameters as 
figure except that the common amplitude is now G = 
50. In the left plot a second atom ò has been coupled to 
the third cavity with a constant coupling 173b = G^t = 5, 
and it is seen that all of the photon ends up in cavity 1 
and 2. Note that atom o is coupled ten times as strongly 
to the field as atom b. In the plot to the right, a furthcr 
third atom c is coupled with a constant coupling <?2c — 
G2C = 5 to cavity 2, and all population now ends up in the 
first state, namely the photon is in cavity 1. These plots 
clearly show how a small disturbance to the adiabatic 
Hamiltonian changes the evolution. If G would have been 
made larger, the perturbations could have been made 
smaller. 



IV. STATE PREPARATION USING ADIABATIC 
TRANSFER AND ATÒMIC MEASUREMENTS 

In the previous sections the atom remained more or 
less in its lower state during the whole process and could 
be seen as an ancillary state, which is never very en- 
tanglcd with the field state. Assuming perfect detec- 
tion efficiency, a measurement on the atòmic state in the 
|±)-basis, after the interaction, would give |— ) with unit 
probability. As long as the atòmic state does not get 
entangled with the field states, an atòmic measurement 
would not modify the cavity states. 

By introducing a third atòmic level \q), which does not 
interact with the field, it is possible to create atom-field 
cntanglement. Thus, an atom in the state \q) will pass 
through the cavities without any interaction, which could 
be due to a large detuning or selection rules. The Hamil- 
tonian is correspondingly only modified by the term for 




-10 -5 5 10 -10 -5 5 10 



t (a.u) t (a.u) 

FIG. 9: This figure shows the dynamics of the same 'star' 
configuration as in figure Q but with small perturbations to 
the Hamiltonian. The left plot shows the same evolution as in 
figure0 but now with a second atom b coupled, in its ground 
state, to the cavity 3. The coupling amplitude G between 
the first atom a and the four cavities are now G = 50, but 
all other parameters are the same as the previous example. 
The coupling between atom fe and cavity 3 is constant during 
the process, g^b = Gíb = 5, and the corresponding detuning 
is zero. The added atom-cavity interaction clearly modifies 
the evolution so that the photon ends up in cavity 1 and 2 
(dashed and dotted lines). To the right we have added yet 
a third atom c, also with a constant coupling Q2 C = C?2c = 5 
and zero detuning, interacting with cavity 2, and now the 
population in that state is removed, so that only cavity 1 is 
populated (dotted line). Note that atom fe and c is much 
weaker coupled to the fields than atom a. Time is given in 
units of T and pulse heights in units of T _1 . 

the atòmic energy in state \q), which could, of course, be 
omitted in a rotating frame. 

In this section we will look at the "H" configuration, 
but other setups could also be considered. We will show 
how it is pos sible to create entangled Schròdinger cat 
states |23l 124 EEj by measuring the atòmic state after 
the interaction. We introduce the atòmic states 

lx)±' b = ^(h> a ' fc ±|ç} a ' b ), (36) 

where the indices a and b refer to the different àtoms. 
We will first couple cavities 1 and 2. From the STIRAP 
evolution 

\0,a,-) — ► |-a,0,-) 

(37) 

|0,Q!,g) — ► \0,a,q) 

for coherent states, it follows, starting from one of the 
atòmic states (|36[1 in the "H" configuration, that 

\0,a,0)\ X ) a + — > i(|-a,0,0)|-) o + |0,at ) 0)|g)°). 

(38) 

After the interaction, the atom is measured in the |x)±~ 
basis, and depending on the measurement result the field 
will be in the state 

JV[|-a,0,0) + (-ir|0,a,0)] , (39) 

where i — for the measurement outcome |x)i an d * = 1 
for the result |x)-j an d the normalisation constant is 
given by N~ 2 = 2[1 + exp(-|a| 2 )]. 
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The atòmic measurement in the desired basis can be 
effected by first using Raman pulses to couple the atòmic 
states |— ) and \q). The resulting unitary evolution should 
transform |x)+ hito |— ) and |x)- hito \q), so that the 
measurement can then be implemented by testing for 
population in the lcvcls |— ) and \q) with a fluorescence 
measurement. With this procedure it is possible to reach 
a very high measurement emciency, almost 100%. Sim- 
ilar methods can be used to implement also generalised 
quantum measurements on àtoms or ions |26j . 

A second atom is then injected into cavity 1 and 3 in 
the state |x)+- The state will evolve into 



^[|0,0, a )|-) b + |- a ,0,0)k) h 
+ (-!)% a,0)\-) b + (-l)%a,0)\q) b ] 



(40) 



Atom b is then measured in the same basis as that for 
atom a, with the result proportional to 

|0,0,a) + (-l) J |-a,0,0) + (-l) ï |0,a,0) + (-l)^'|0,«,0), 

( 41 ) 

for the cavity field states, where j is defined as i is, but 
for atom b. We have here left out the normalising con- 
stant, since it will depend on the measurement outcome 
for atom b. Depending on the known measurement out- 
comes for àtoms a and 6, we are able to prepare four 
possible entangled states, 

|*oo) oc (| - a, 0, 0) + 2|0, a, 0) + |0, 0, a)) , i = j = 

|*oi) cx (-| - o, 0, 0) + |0, 0, a)) , i = 0, j = 1 

|*io) cx (| -a,0,0)-2|0,a,0) + |0,0,a)), i = 1, j = 

|*n) cx (_| _ a ,0,0) + |0,0,a)) , i = j = 1. 

(42) 

We may also consider the following scenario. If the 
second atom is injected in the state \—) b instead, it will 
leave the setup in the same state, and the resulting field 
state is 



N(\0,0,a) + (-l)*|0,a,0)) 



(43) 



Let us fix (3 through a = 2(3 and introduce the displace- 
ment operator D with the properties 



D(fi)\a) = e dm ^^\a + /3). 



(44) 



If the operator D(—/3) is applied to both cavity 2 and 
3 and for real a and f3, the resulting entangled state of 
cavities 2 and 3 becomes 



N (\ — (3, 0) + (— 1)'|/3, —(3)) , 



(45) 



where N is defined as before. Here cavities 2 and 3 are 
both in a Schròdinger cat state and entangled with each 
other. This kind of entangled state is of great interest for 
quantum teleportation |27| and quantum computing with 



coherent states |28|, but also for studying quantum phe- 
nomena in general, likc cntanglement and decoherence in 
the classical limit |2íj. Using a 50/50 beam splitter, this 
state may be transformed into \\f2(3, 0) + (— 1) ! | — V%/3, 0), 
i.e. a cat state in one of the modes only, with vacuum in 
the other mode. 

It should be mentioned that the atòmic states \x)± 
could have been defined in different ways, leading to other 
entangled field states. The initially prepared and mea- 
sured atòmic basis need not be the same. We could have 
considered different setups of cavities and àtoms and the 
initial coherent state could have been any state, for ex- 
ample squeezed states. 

We conclude this section by considering another exam- 
ple of how to prepare a Schròdinger cat state. We now 
assume just two overlapping cavities and a single atom 
as in section [H] The difference is that the atom a now 
should have (at least) two degenerate ground state levels 
|— )i.n labeled by I and II, such that the coupling ampli- 
tudes are G la j = G 2a ,i = Gi a ,u = G and G20.11 = -G, 
where 1 and 2 indicate the cavity and 1,11 the transition. 

One way to achieve this might be to impose a chosen 
quantization axis for the atom using an external elèctric 
field, thus forcing the dipole moment d of the atom to 
have the suitable components along the directions of the 
two làser fields. Alternatively it may be possible to use 
sclcction rules for the transitions in such a way, that it 
is possible to choose the signs of the elèctric field com- 
ponents inducing the different transitions. The choicc 
should be made in such a way that d • E has the re- 
quired signs for the four different combinations of làser 
and atòmic transition. 

Assuming that this choice of coupling constants is pos- 
sible, if we now prepare the atom in state |— )f , an initial 
coherent state \a, 0) in mode 1 will be transferred into 
|0, — a) in mode 2. This is because as we can see from 
the discussion in section Hi Al when Gi a j/G2a,i > 0, then 
an arbitrary field state /(àj)|0) in cavity 1, will be trans- 
ferred into a state /(— a|)|0) m cavity 2. But if the atom 



is prepared in 



/in 



an initial coherent state \a, 0) in 



mode 1 will be transferred into |0, a) in mode 2, without 
the minus sign. Again, this is because Gi 0j n/G2 a ,ii < 0, 
so that an arbitrary field state in cavity 1, f(à\)\0), will 
be transferred into a state /(ò^lO) in cavity 2. If the 
atom is initially in a superposition of the two states, 
= l/V2(|-)f ± the result will be 

1 



V2 



(|0, -«>|->f ± |0, «>!->&) . (46) 



This is a Schròdinger cat state for cavity 2 and the atom. 
If we wish to disentangle the atom and the cavity, the 
atom may be measured in the basis l/y/2(\— )" ± |— 
Depending on the measurement outcome, we are left with 
one of the states 



N(\0,a) ± |0,-o)). 



(47) 



The coherent state is transfered from cavity 1 into a cat 
state in cavity 2. 
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V. CONCLUSIONS 

In this paper, we have given several examples of cav- 
ity field state preparation and transfer using adiabatic 
methods. The technique we use is related to stimulated 
Raman adiabatic passage (STIRAP) 0,0. In Standard 
STIRAP, atòmic energy levels are coupled by làser pulses 
in order to transfer population between the atòmic states. 
In the present scheme, cavity field mode are effectively 
coupled by àtoms in order to transfer population between 
the cavity modes. A previous paper showed that not only 
photon number states, but arbitrary cavity field states 
can be transferred using this method In this paper, 
we have in particular considered preparation of entan- 
gled states of two or more cavities, such as an EPR state 
and a W state, and various entangled superpositions of 
coherent states in different cavities. The theoretical con- 
siderations are supported by numerical simulations. It 
may also be possible to use similar techniques in sòlid 
state systems, replacing the cavities and àtoms in our 
discussion with cavities coupled to Josephson junctions 

Eü). 

One advantage of adiabatic state transfer and prepa- 
ration methods is that they are relatively robust against 
changes in the individual coupling pulsc strengths and 
pulse durations. In contrast, state transfer e.g. in the 
Jaynes-Cummings model [Ï5| relies on the ability to ex- 
perimentally control the areas of coupling pulses very ac- 
curately. The situations considered in this paper are by 
no means totally unrealistic considering the present sta- 
tus of experiments in QED. An important condition is 
that all the cavity modes have to be degenerate. This 
results from energy conservation; if the modes were not 
degenerate, perfect state transfer between modes would 
not be possible. The adiabaticity for processes like the 
ones considered in this paper, is roughly given by the 



coupling amplitude times the pulse width Ger, see 1]. In 
the example of figure |21 we have Ga — 300, while us- 
ing typical experimental vàlues of G/2ir ~ 100 MHz and 
a ~ 0.3 s _1 H,[^ï|, the adiabaticity parameter becomes 
Ga w 200. With these characteristic non-scaled parame- 
ters, the coupling is multiplied by 27r • 10 6 and the time 
seales by 10" 7 , the adiabatic transfer of figure [3 gi ves 
a final population of 96.9 % in the target cavity, while 
the maximum population in the intermediate cavity 1, 
during the process, is 2.2 %. The whole operation, with 
the two àtoms passing through the three cavities, takes 
about 2 /js, which is much shorter than the characteristic 
life times of cavity and the atòmic states [U- Remem- 
ber, that the example of adiabatic transfer in figures [21 
and |21 involves more virtual intermediate levels than, for 
example in the generation of EPR and W states in fig- 
ures and Another point to emphasize is that, as 
the field amplitude increases, the number of intermedi- 
ate states also increases, which makes the adiabaticity 
constrains stricter. Even though it is possible to have 
strong enough couplings and small decay rates for realiz- 
ing the schemes proposed in this paper, it is not obvious 
whether it is a simple task to add crossing cavities in 
current experimental setups. 

In conclusion, adiabatic techniques offer rich possibili- 
ties for state transfer and population. 
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